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Abstract. The three-dimensional strongly screened vortex-glass model is studied 
numerically using methods from combinatorial optimization. We focus on the effect 
of disorder strength on the ground state and found the existence of a disorder- 
driven normal-to-superconducting phase transition. The transition turns out to be 
a geometrical phase transition with percolating vortex loops in the ground state 
configuration. We determine the critical exponents and provide evidence for a new 
universality class of correlated percolation. 
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1. Introduction 

The gauge glass model is a paradigmatic model for disordered arrays of Josephson 
junctions or amorphous granular superconductors It has been argued that it also 
describes the relevant physics of the superconductor-to-normal phase transition in high- 
Tc superconductors p|. A powerful tool to investigate this transition is the domain 
wall renormalization group (DWRG) technique that has been applied successfully to 
this model P, H: In essence one calculates the stiffness of the system with 

respect to twisting the phase variables at opposite boundaries of a system of linear 
size L. If the twist costs an energy that increases with L one concludes that the system 
is superconducting, if it decreases, one concludes that phase coherence necessary for 
superconductivity is destroyed by thermal fluctuations, i.e. the system is in a normal 
phase. In this paper we study this model in the strong screening limit with varying 
strength of the disorder at zero temperature. We will find a superconductor-to-normal 
transition (at T = 0) at a critical disorder strength and show that it is accompanied by 
a proliferation of disorder induced global vortex loop. By a finite-size scaling analysis 
of the loop statistics we show that it is a percolation transition of a novel universality 
class. 

This paper is structured as follows: in section |^ we present the model a motivation 
to expect a disorder-driven phase transition using the concept of the defect energy. In 
the next both sections our results are presented. Section 3 shows a clear phase transition 
via the study of an excitation loop perturbation. In section 4 the transition is shown 
to be a geometrical phase transition, what gives rise to apply percolation theory to the 
vortex glass model. The critical probability, above which a loop percolates, the critical 
exponents and scaling relations are calculated numerically. We close with a summary 
in section ^ 

2. Model 

The phenomenological lattice model describing the phase fluctuations in a granular 
disordered superconductor close to the normal-to-superconducting phase transition is 
the gauge glass model 0, |^ 

n = -Jj] cos(0i - 0, - Aij - \-^aij) + \ ^(V X a)2, (1) 

where J is the effective coupling (set to 1) and (pi the phase on site i. The sum is over all 
nearest neighbors {ij) on a simple cubic lattice of system size L with periodic boundary 
conditions. Aij are the vector potentials, which are uniformly distributed on 

Aij e [0,2na] with a fixed a E [0, 1], (2) 

where a defines the disorder strength, a = 1 corresponds to strong disorder and a = 
to the pure system, respectively. A is the bare screening length. The fluctuating vector 
potentials a^j are integrated over from — oo to oo subject to V • a = 0. The last term in 
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(|^) describes the magnetic energy and its sum is over all elementary plaquettes of the 
lattice. To investigate the gauge glass model in the strong screening limit A ^ we 
make use of the vortex representation 0, which gives after standard manipulations 

Hy^^ = - ~ with the magnetic field bj = — ^ Aij (3) 

2 . ^TT g 

subject to the local constraint (V ■ n)j = 0. The computation of the ground state 
of the Hamiltonian (^), i.e. the vortex configuration n with the lowest energy 7Yy(iij, 
is a minimum-cost-fiow problem that can be solved exactly in polynomial time with 
appropriate combinatorial optimization algorithms |^. 

We use the defect energy scaling method to show that there is a superconducting- 
to-normal phase transition at low temperature T varying the strength of disorder a. The 
idea is to calculate the energy AE necessary to introduce a low-energy excitation loop 
(or domain wall) of size L to the system. We generate the excitation loop by a global 
manipulation of the energy couplings along a fixed direction, as described in reference 
0, 3 in detail. The defect energy AE results from the difference energy of ground state 
with and without a global excitation loop. Its disorder average is assumed to scale with 
the system size L as 

AE ~ (4) 

The sign of the stiffness exponent 6 determines whether the ground state is stable 
with respect to thermal fluctuations. If > it costs an infinite amount of energy to 
induce a domain wall crossing an infinite system {L oo) and therefore the the ground 
state remains stable at small but non-vanishing temperatures: there is an ordered low 
temperature phase, like in a 3d XY-ferromagnet. On the other hand, if 6' < arbitrarily 
large excitations loops cost less and less energy: the ground state is unstable and thus 
not an ordered phase at any non-vanishing temperatures, like in a 2d XY spin glass. 

We can easily see from (§) that for small disorder (i.e. small a) the ground state is 
simply n = and 6=1: For a given disorder strength a, it is bi G [—2 a, 2 a]. Thus, as 
long as o" < 1/4 it is 1 6j I < 1/2 and the absolute minimum of all terms (11^ — 64)^ occurring 
in the Hamiltonian (^) fulfilling the constraint that rij has to be integer is rij = 0. A 
global excitation loop contains at least L bonds with rij = 1, which implies that for 
(T < 1/4 an additional global excitation loop would cost a defect energy AE oc L. This 
implies a stiffness exponent ^ = 1 for small disorder, certainly for a < 1/4, possibly even 
for larger sigma, as we will see below. Thus we can already at this point conclude that 
the for weak disorder the system described by (^ is superconducting (or ferromagnetic 
in the the magnetic, XY language), as it is in the pure case (a = 0). 

On the other hand, in the opposite limit of strong disorder, a = 1, defect energy 
calculations |^ gave a negative stiffness exponent 9 = —0.96 ± 0.05, indicating the 
absence of an ordered low-temperature phase (in particular the absence of a stable low 
temperature vortex glass phase [rn||). Therefore one can expect a disorder driven phase 
transition at zero temperature from a superconducting phase for weak disorder to a 
normal phase for strong disorder. We expect that this transition takes place at a critical 




Figure 1. Log-log plot of the disorder averaged defect energy /S.E vs. system size L 
for different disorder strengths a. For L = 4, 6, 8 we used Ngamp = 20000 samples, for 
L = 12 5000 samples, L — 16 1000 samples and L = 24 500 samples, respectively. 



disorder strength dc (cTc > 1/4 from what we said above) and is characterized by a 
discontinuous jump of the stiffness exponent 6 from 1 to —0.96 (here we assume the 
simplest scenario in which one has only two attracting zero temperature fixed points 
besides the critical point (Xc)- 

3. Defect Energy 

Figure |I], showing the defect energy AE versus system size L in a log- log plot, 
demonstrates that our numerical results confirm this hypothesis. The slopes of the 
different curves, representing different disorder strengths a, are identical to the stiffness 
exponent 6. We observe that around cTc = 0.495 ± 0.005 it jumps from positive to 
negative with increasing a, this is our estimate for the location of the disorder driven 
transition from the superconducting to the normal phase. Note that for the unscreened 
gauge glass XY model it was found (7^ ~ 0.55 

In what follows we will show that this zero temperature transition is actually a 
2nd order phase transition characterized by a single length scale diverging at cTc. This 
length scale corresponds to the average diameter of closed loops in the ground state and 
at (Tc these loops percolate the infinite system (note that we have periodic boundary 
conditions), see Fig. |[ Thus what actually happens at Uc is a percolation transition of 
vortex loops in the ground state and in order to estimate its critical exponents we will 
perform a finite-size scaling analysis of the critical behavior now. 
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Figure 2. Three ground state configurations for growing disorder strengths a = 0.45 
(left), 0.50 (middle) and 0.55 (right), respectively, and system size L = 16. At ct = 0.50 
it appears a percolating loop (dotted line). 



4. Vortex Loop Percolation Transition 

First we note that at cTc the concentration of vortex variables that are non-zero (n^ ^ 0), 
i.e. the probability p with which a bond in the simple cubic lattice is occupied with 
a vortex segment, turns out to be Pc = 0.033 ± 0.005. This value is much lower 
than the percolation threshold for conventional bond percolation on the simple cubic 
lattice |11|], which is piP'^^'^'^ ^ 0.249 [jl2|. This is a consequence of the global constraint 



(divergence-free) underlying the optimization problem (^, which obviously causes strong 
correlations in the bond occupation process. Hence we suspect that the transition we 
are considering establishes a new percolation universality class p!3| . 

The geometrical objects of the ground state n of model (^ that we are going to 
study are loops. The algorithm to detect loops is the following: 
given the ground state configuration n; 

while it exists a vortex segment with n^ 7^ along the bond i do: 

(I) choose i] 

(II) find the shortest path Pj along non-zero vortex segments from the target site of i 
to the source site of i, where the direct path along i is excluded; 

(III) calculate flow := Z]jePiU{i} '^j- if flow 7^ or if there are two occupied bonds in a 
distance L or larger along the x, y or z direction, which belong to the same loop, then 
the loop is called a global loop else a local loop; 

(IV) cancel the detected loop Pi U {i} and continue the while loop. 

For each system with L = 4, 6, 8, 10, 12, 16 we calculated 2000 different samples, 
for L = 20, 24 and 32 1000 samples, respectively, and then analyzed the loop statistics. 

By studying the probability Pi^^'"^°{a) that a system of linear size L contains at 
least one percolating loop we can check, whether the percolation transition does indeed 
coincides with the jump in the stiffness exponent located above. Its finite-size scaling 
form is given by 

Pf ^^°(p) = pP--[ (a - a,) ■ L^/" ] , (5) 
thus it is system size independent at Oc and curves for different system sizes should 
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Figure 3. Finite-size scaling of the percolation probability pp^''=° (left) and the 
average number Np^rco of percolating loops (right) for different system sizes L with 
CTc = 0.492 and v — 1.05. The inset shows the raw data. The error bars are smaller 
than the symbol size and therefore omitted. 



intersect. Our data are shown in the inset of Fig. ^ (left) and we locate the intersection 
point 

= 0.492 ± 0.005 (6) 

agreeing well with our estimate for from the defect energy analysis. 

Next we deduce an estimate for the correlation length exponent v by plotting 
ppercoi^^-^ versus (cj — (Tc) ■ L^^'^ ^ where we fix cXc and determine v such as to achieve 
the best data collapse. This is done in Fig. |^ (left) and we obtain 

1/ = 1.05 ±0.05. (7) 

This estimate for v lays between the value of the two- and three-dimensional bond 
percolation [|lT] . 



The analysis of the average number of percolating loops N^~^'^°[p), obeying a similar 
finite-size scaling form N'j^^'^°{p) = Np^^^"^ (cr — cr^) . L^/'^ ] gives the same estimates for 
(Jc and z/, c.f. Fig. |^ (right). Note that at the average number of percolating loops 
does not (or only weakly) depend on the system size and is small: Nf^~^'^°{pc) ~ 0.3. The 
maximum number of percolating loops we observed for L = 32 at was 3 with a very 
low probability. 

The average mass m of a percolating loop at (Jc scales with L like, 

m~L^^ (8) 
where is a fractal dimension. For a = cXc we get with the data shown in figure ^ 

c/j = 1.64 ±0.02 . (9) 

The probability Pqo that a bond belongs to a percolating loop is expected to scale 

hke 

P^r^L-P/'' P^[{a-a;)L^I% (10) 



Superconductor-to- Normal Phase Transition in Vortex Glass Model 



7 




Figure 4. Left: Plot of the average mass m of a percolating loop vs. L ioi a — Uc- 
The error bars are smaller than the symbols. The straight line is a least square fit to 
m ^ ^ giving df = 1.64±0.02 Right: Finite-size scaling plot of the probability Pqo 
for a bond belonging to a percolating loop for different system size L with tXc = 0.495, 
1/ — 1.05 and (3/v — 1.4. The inset shows lin-log plot of the raw data. 
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Figure 5. Left: Probability distribution Um at ct = CTc for different system size L. 
A least square fit to n„i ~ m~r yields r = 2.8 ± 0.1. Right: Finite-size scaling of rim 
for L = 32 with tJc = 0.495, s = 0.6 and r = 2.95. For m > 30 the statistics is over 
less than 1000 loops for each a. 



The figure ^ (right) shows the raw data of Poo (inset) and the plot of the scaling law 
(|ig) with u = 1.05 ± 0.05 and p/u = 1.4 ± 0.1, i.e. 

/3=1.4±0.1. (11) 

The usual hyper-scaling relation, jS/u = d — df, known from conventional percolation 
111 gives df = 1.6±0.1, which is consistent with (^. 

The loop distribution function n^, i.e. the average number of finite loops of 
mass m per lattice bond obeys the scaling form (in the limit L oo) 

Um ^ m'^ Umiio- - ac)m''), (12) 
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Figure 6. Left: Finite-size scaling plot of the susceptibility x for = 0.495, 
V ^ 1.05 and ^/v = 0.42 Right: Scaling plot of the correlation length ^ for CTc = 0.495 
and v — 1.05. The inset shows the raw data. 



where r is the Fisher exponent and s another critical exponent (usually denoted a 
in conventional percolation, which we avoid due to possible confusion with the disorder 
strength a). The exponent s describes how fast the number of loops of mass m decreases 
as function of m close to Uc. Figure ^ (left) shows the raw data of for different L and 
a = cTe. For L = 32 we get r = 2.89 ± 0.05 and for L = 127 (3 samples) r = 2.84 ± 0.06, 
respectively. In the limit L — > oo we expect 

r = 2.8 ±0.1. (13) 



From the finite-size scaling plot of equation (|T2D, we get r = 2.95 ± 0.05 and 

5 = 0.6 ±0.1 (14) 

for = 0.495 ± 0.005 and L = 32 in figure | (right). 

The zeroth moment n = J2m represents the average number of loops per bond. 
Below (Jc the data collapse and satisfy n ^ a. The average loops size, defined as the 
ratio of the second and first moment of the loop distribution [|11|]: 

X-=iYl T^^^rn \ I mUm] ■ (15) 

\m=4 / \m=4 / 

is expected scale like 

X-^L^''' x[{a-a,)-L'/']. (16) 

The data in figure |] (left) show the raw data (inset) and verify the scaling law ([T6|) with 
7/z/ = 0.4 ± 0.1 for ac = 0.495 and u = 1.05 ± 0.05, i.e 

7 = 0.4 ±0.1. (17) 

The above estimates for 7 ([T7|) and (3 ([TT|) together with those for r (|13D and s (0) 
fulfill the usual exponent relation known from conventional percolation 

, = . ;3 = ^. (18) 
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Figure 7. Left: Average radius of gyration Rm vs. the mass m ai a = cjc for 1000 
samples and L = 32. A least square fit to Rm ^ fn'' yields p — 0.66 ± 0.02. Right: 
Finite-size scaling of the square loop radius of for different system sizes L with 
CTc — 0.492, 1/ = 1.05 and xjv — 0.8. The inset shows the raw data. 



Near Oc the linear size of a finite loop is characterized by the correlation length ^, 
which we calculate with the help of the radius Rm^ of gyration for the loop i of mass nii 
defined as 

•= — E I'^ji - J^o^r with ro, := — E ""Ji ' (19) 

iil'i j—i iil'i j—i 

where rji is the position of a bond j of the loop i and roi the center of mass, respectively. 
Then, the correlation length ^ is defined by 

:= f E RWnm) I f E ^'^™) ' (20) 

\m=4 / \m,=4 / 

where Rm is the average radius of gyration of loops of mass m (averaged over disorder 
and individual loops). The raw data of ^ are shown in inset of figure |^ (right). The 
finite-size scaling form for ^ is 

i^L.l\{a~a:).L'l^\. (21) 

From the best data collapse we get v = 1.05±0.05, as shown in figure ^ (right), consistent 
with (g). 

At the percolation threshold the average radius of gyration Rm of a loop of mass 
m increases algebraically 

Rm ~ mP. (22) 

In figure |^ (left) we plot Rm for a = ac and fit the data in the interval m G {10, 100} 
to the power law (|22|) , which yields 

p = 0.66 ±0.02 or = 1/p = 1.51 ± 0.05 , (23) 

which agrees with our previous estimate for the fractal dimension dj oi the percolating 
loops (^ within the error bars. 
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Another quantity, which characterizes the size of finite loops, is the mean square 
radius i?^, defined as 

:= ( E C ^ ) / ( E m UrA . (24) 

\m=4 / \m=4 / 

We expect to scale like 

i?2_L"/"i?[(a-ae)Li/n> (25) 

where x is another critical exponent. As depicted in figure |^ (right) for the best data 
collapse we get x/v = 0.8 ± 0.1 with u = 1.05 ± 0.05 (and cxc = 0.495 ± 0.005), i.e. 

x = 0.8 ±0.1. (26) 



This exponent should fulfill the relation [11 



x = 2iy-P (27) 
With u from (|^ and f3 from (|Tl|) we get x = 0.7 ± 0.2, which is consistent with (|26|). 



5. Summary 

In summary, we studied the ground state of the three-dimensional strongly screened 
vortex glass model, numerically. We found a clear evidence for a disorder- driven 
superconducting-to-normal phase transition indicated by a change in the stiffness 
exponent at cTc- This transition turned out to be a percolation transition for disorder 
induced vortex loops crossing the whole system. 

At first sight it might be surprising, why the existence of percolating vortex loops 
is related to a change in the stiffness exponent of model (|I]). However, the the stiffness 
exponent provides information on how hard it is to induce a domain wall into a system 
of linear size L and a domain wall is surrounded by a global vortex loop. If, at and 
above a critical disorder strength, global vortex loops proliferate already in the ground 
state, the creation of an extra excitation loop will, with probability one, costs only an 
infinitesimal amount of energy in the infinite system size limit. 

A similar observation — the coincidence of vortex loop percolation and a thermal 
phase transition in superconductors — has been made earlier in models for high-Tc 
superconductors: In [|I^ it was shown for a model of a pure superconductor that the 
melting transition of the Abrikosov fiux line lattice at the temperature Tc2, where the 
transition from the superconductor to normal phase takes place, is accompanied by a 
proliferation of thermally induced global vortex loops. And similarly in |jl5| it was shown 
that the temperature driven resistivity transition in disordered high-Tc superconductors 
is also accompanied by a percolation transition of vortex lines perpendicular to the 
applied field. These thermally induced transitions are, however, in universality classes 
different from the disorder induced transition we studied here. 
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